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REMARKS 

The Office Action of February 17, 2006 has been received and its contents carefully 
considered. 

Claims 1 to 5 are all the claims pending in the application. 

The Examiner advises that he has interpreted the phrase "3d transition metal elements" 
used throughout the claims to be limited to the elements Cr, Mn, Fe, Co, Ni and Cu as set forth in 
Paragraph [20] of the specification. 

In response, applicants have amended claim 1 to recite the metal elements described in 
Paragraph [20] of the application. 

Claim 1 has been rejected under 35 U.S.C. § 102(b) as anticipated by EP 0 364 631 Al to 
Mizimioto et al. 

Applicants submit that EP '631 does not disclose or render obvious the presently claimed 
invention and, accordingly, request withdrawal of this rejection. 

The present invention as set forth in claim 1 as amended above is directed to a magnetic 
crystalline alloy comprising Pt in an amount of 40 at% to 60 at%, and at least two 3d transition 
metal elements selected from the group consisting of Cr, Mn, Fe, Co, Ni, and Cu, wherein the 
total amount of the 3d transition metal elements is from 60 at% to 40 at%, and the average 
number of valence electrons in the respective 3d transition metal elements as calculated on the 
basis of the compositional proportions of the elements is from 7.5 to 8.55. 
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Thus, applicants have amended claim 1 to recite that the alloy is crystalline, and that the 
average number of valence electrons in the respective 3d transition metal elements as calculated 
on the basis of the compositional proportions of the elements is from 7.5 to 8.55. 

Support for the amendment to claim 1 to recite the value 8.55 can be found in the present 
specification at page 13, Table 2, Example 2. 

In the Response Under 37 C.F.R, § 1.111 filed on December 2, 2005, applicants pointed 
out that the average number of valence electrons in the respective 3d transition metal elements of 
the Pt46NiiCo53 alloy in Example 3 of EP *63 1, upon which the Examiner relied, would be 
calculated to be 9.02. 

In the present Office Action, the Examiner states that since the calculated 9.02 value is 
very similar and closely approximates the claimed upper limit of 9 of original claim 1 , one of 
ordinary skill in the art would have expected the example alloy of EP *631 and the presently 
claimed alloys to have the same properties. 

In response, and as discussed above, applicants have amended claim 1 to recite that the 
upper limit is 8.55. 

Since claim 1 recites that the upper limit be 8.55 and the alloy the Examiner has 
identified has a value of 9.02, appUcants submit that the Pt46NiiCo53 of EP '63 1 does not 
anticipate claim 1. 

Further, applicants again point out that technology field is different between the present 
invention and EP '63 1 . As described in the "Object of the Invention" on page 3 of EP '63 1 , the 
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invention of EP '631 is basically related to the technology field of a magneto-optical recording 
media, recently often called "MO", which is totally different from that of the present invention. 

The present invention relates to a magnetic alloy for a magnetic layer used in a magnetic 
recording medium for a so-called "hard disk" device. Accordingly, the characteristics necessary 
to the present invention are totally different from those of the magneto-optical recording media, 
and do not need the parameters disclosed in Table 1 on page 9 of EP '63 1 . 

In addition, EP '63 1 discloses that an amorphous alloy is necessary for a magneto- 
optical recording media. See the entire disclosure of EP '631, including the Title, Field of 
Invention Object of the Invention, and Summary of the Invention, which indicate that the 
invention of JP '63 1 is directed to an amorphous alloy. However, the present invention relates to 
a magnetic alloy having a crystallized structure, which is obvious to a person having ordinary 
skill in the art who reads paragraph [25] at page 7 of the present specification describing the 
order parameter, or the orientation of the crystal such as (001) or (002). 

Accordingly, applicants submit that the reasons for the rejection are xmfounded. 

In particular, applicants submit that the Examiner's statement that Example 3 at page 9 of 
EP '631, which is an amorphous alloy, anticipates claim 1 of the present application is in error. 

With respect to applicants' arguments that the EP '631 alloy is amorphous whereas the 
alloy of the present invention has a crystalline structure, the Examiner states that he is not 
persuaded because the present claims are silent with respect to the crystalline structure of the 
claimed alloy. The Examiner states that the present claims, therefore, encompass the amorphous 
structure taught by EP '63 1 . 
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In response, applicants have amended claim 1 to state that the alloys of the present 
invention are crystalline. Support for this amendment to recite a crystalline structure for the 
claimed alloy is found in the present specification in its description of the alloy. 

For example, paragraph [22] of the present specification states that if the total amount of 
the 3d transition metal elements exceed 60 at%, "the structure of the magnetic alloy changes Llo 
to another structure". A "Llo" structure refers to a particular crystal structure among the many 
crystal structures that are known to exist. See the attached Attachment A, which is a compilation 
of various printouts that were obtained from the following website links: http://cst- 
www.nrl.naw.mil/lattice/index.html ; http://cst-www.nrl.navv.mil/lattice/struk/index.html ; 
http://cst'Www.nrl.navv.mil/lattice/struk/ltvpe.html ; and http://cst« 

www.nrl.naw.mil/lattice/struk/l 1 O.html . These printouts describe "Crystal Lattice Structures." 

As can be seen from the enclosed Attachment A, crystal lattice structures can be 
designated in a variety of different ways, including a "Strukturbericht Designation." As set forth 
in Attachment A, Strukturbericht symbols (designations) are a partly systematic method for 
specifying the structure of a crystal. The L designation is used to describe alloy crystal types, 
and among the Strukturbericht Type L structures, there is the Llo (AuCu) structure. 

Thus, the reference to the Llo structure in the present specification clearly shows that the 
alloys of the present invention are crystalline. 

Moreover, in the claims, (001) and (002) refer to a "plane direction", which is used to 
describe a crystal structure, and which is never used to describe an amorphous structure. See, for 
example, Fig. 2 on page 1707 of the attached Attachment B, Le Bouar, Y., "MECHANISM 
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AND MODELING OF SAW-TOOTH STRUCTURE FORMATION IN THE Lh-Llo TWO- 
PHASE SYSTEM," Acta mater.. Vol. 48 (2000), pages 1705-1717. 

EP '631 shows in Example 3 of Table 1 that the Pt46NiiCo53 alloy is amorphous for its 
crystalline condition, and its Kerr rotation angle (0kO) is 0.15 and Hco (coercive force) is 0.9 for 
the magnetooptical characteristics. EP '631 also shows in Example 22 that the alloy Pt75Co25 is 
crystalline for its crystalline condition, but its Kerr rotation angle (6k0) is unmeasurable, and 
Hco (coercive force) is also unmeasurable . This means that the characteristics or purpose of the 
amorphous alloy of EP '631 cannot be applied to a crystal alloy. 

On the other hand, the magnetic alloy of the present invention has a crystal structure, and 
does not have the same characteristics as the magnetooptical effects of the alloy of EP '631, such 
as Kerr rotation angle (GkO). 

In view of the above, applicants submit that EP '63 1 does not anticipate claim 1 and, 
accordingly, request withdrawal of this rejection. 

Claims 2 and 3 have been rejected under 35 U.S.C. § 102(b) as anticipated by or, in the 
alternative, under 35 U.S.C. § 103(a) as obvious over EP '631. 

The Examiner states that EP '631 is silent with respect to the properties recited in claims 
2 and 3. The Examiner argues that the alloy disclosed in EP *63 1 would be expected to possess 
all of the same properties as recited in present claims 2 and 3 because the alloy disclosed in EP 
'63 1 has a composition that is encompassed by the present claims. 

In response, and as discussed above, applicants submit that the alloy of Example 3 of EP 
'63 1 would not be expected to possess all of the same properties as recited in claims 2 and 3. 
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As stated above, the characteristics necessary to the magnetic alloy of the present 
invention are totally different from that of the magneto-optical recording media of the EP '631, 
and the magnetic alloy of the present invention do not need the parameters disclosed in the Table 
1 onpage9ofEP'631. 

In the Response Under 37 C.F.R. § 1.1 1 1 filed on December 2, 2005, applicants stated 
with respect to the "S" value recited in claim 2, that "although it may be possible to calculate the 
*S' value for the alloy of the Example 3 of EP '631, such value *S' will be zero, which has no 
technical meaning, because of the characteristics of the amorphous alloy." 

The Examiner states in the present Office Action that applicants have not provided any 
evidence to support the assertion that the "S" value of EP '631 is 0. 

In response, applicants provide the following clarifying remarks. 

In particular, applicants submit since the alloys of EP '631 are amorphous, the "S" value 
of the alloy of EP '631 cannot be obtained and, therefore, it is meaningless to refer to an "S" 
value of an amorphous alloy. Amorphous alloys do not possess an "S" value. This can be 
illustrated by referring to the attached Attachment C containing X-ray diffraction profiles shown 
in Figs. 1-3, which Mr. Kanazawa, one of the inventors, has provided. 

Figs. 1-3 of Attachment C show X-ray diffraction profiles for Fe-Co-Pt alloys. As 
described in the present invention, the order parameter "S" is calculated from the integration 
intensity (1(001)71(002)). 

In order to calculate S, at the least, the X-ray intensity from the (002) plane needs to be 
observed, such as in Fig, 1 and Fig. 2 of Attachment C, which means these alloys have a crystal 
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structure. However, once the structure becomes amorphous, (001) and (002) peaks cannot be 
observed, as is shovm in Fig. 3 of Attachment C. This means that the "S" value cannot be 
defined for an amorphous alloy, and 1(001) cannot be observed. Therefore, the amorphous alloys 
of EP '631 do not possess an "S" value according to claim 2 of the present application. 

Further, it is impossible to calculate the Ku value as set forth in claim 3. 

In view of the above, applicants submit that EP '63 1 does not disclose or suggest the 
subject matter of claims 2 and 3 and, accordingly, request withdrawal of this rejection. 

Claims 4 and 5 have been rejected under 35 U.S.C. § 103(a) as obvious over EP '63 1. 

The Examiner states that EP '631 teaches that the alloy of Example 3 can be used as a 
magnetic memory medium as disclosed at page 8, lines 46-49. 

Applicants continue to believe the Examiner intended to refer to page 6, lines 46-49 of 
EP '63 1 , where it is disclosed that the amorphous alloy thin films have an easy axis of 
magnetization perpendicular to the film face and can be used in such various fields as magnetic 
recording materials, including vertical magnetic recordmg films, magnetic bubble memories and 
magnetooptical recording. 

Claims 4 and 5 are dependent claims that depend ultimately from claim 1 or 2. 

Accordingly, applicants submit that claims 4 and 5 are patentable for the same reasons as 
discussed above with respect to claims 1 and 2. 

As stated above, the material itself of the present invention having a crystal structure is 
totally different from that of EP '631 having an amorphous structure. 

In view of the above, applicants request withdrawal of this rejection. 
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In view of the above, reconsideration and allowance of this application are now believed 
to be in order, and such actions are hereby solicited. If any points remain in issue which the 
Examiner feels may be best resolved through a personal or telephone interview, the Examiner is 
kindly requested to contact the imdersigned at the telephone number listed below. 

The USPTO is directed and authorized to charge all required fees, except for the Issue 
Fee and the Publication Fee, to Deposit Account No. 19-4880. Please also credit any 
overpayments to said Deposit Account. 

Respectfully submitted. 



SUGHRUE MION, PLLC Sheldon I. Landsman 

Telephone: (202)293-7060 Registration No. 25,430 

Facsimile: (202) 293-7860 

WASHINGTON OFFICE 

23373 

CUSTOMER NUMBER 

Date: May 17,2006 
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Crystal 
Lattice-Structures 
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Strukturbericht Designation 


Pearson Symbol 


Space Group 


Prototype 


Structures of Intermetallic Alloy Phases 



This page offers a concise index of common crystal lattice structures. A graphical representation as 
well as useful information about the lattices can be obtained by clicking on the desired structure 
below. 

This page currently contains links to 254 structures in 90 of the 230 space groups. 
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Sulfur and its Compounds 
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and related structures 
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and Related Structures 
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and Related Structures 
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Other Structures 



Structures indexed by: 

• Strukturbericht Designation 

• Pearson Symbol 
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• Prototype 



References used in creating these pages. 
Programs used in creating these pages. 



Other sources of information about crystal lattices can be found at various sites around the Web . 



NRL policy prohibits mention of the creators of this page. 



Questions, comments, additions, and corrections are welcome. r Privacy Advisory ) but first please 
read the FAQ. If you would like to see a structure added to the database, please read here first . 

Return to the Code 6390 Home Page. 



This page maintained by lattice@dave.nrl . navv.mil (Privacy Advisory) 

Last updated: 21 October 2004. 

Current URL: http;/ /est -www. nrl.naw.mil /lattice/index. html 

The appearance of external hyperlinks does not constitute endorsement by the United States Department of Defense, the United States Department of the Navy, and the Naval Research 



http://cst-www.nrl.navy.mil/lattice/index.html 



5/16/2006 



Crystal Lattice Structures 



Page 4 of 4 



Laboratory of the linked web sites, or the information, products or serveices contained therein. For other than authorized activities such as military exchanges and Morale, Welfare, and 
Recreation (MWR) sites, the United States Department of Defense, the Department of the Navy, and the Naval Research Laboratory does not exercise any editorial control over the 
information you may fmd at these locations. Such links are provided consistent with the stated purpose of this DoD web site. 



Return to CCMS home page 



http://cst-www.nrl.navy.mil/lattice/index.html 



5/16/2006 



Crystal Lattice Structures: Index by Strukturbericht Designation Page 1 of 2 

BEST AVAILABLE COPY 



Sites:-' 

m 



Crystal Lattice Structures: I^^JiZ'Z 

Index by Strukturbericht Designation 



• Strukturbericht symbols are a partly systematic method for specifying the structure of a 
crystal. Thus the A structures are supposed to be monatomic, B's are diatomic with equal 
numbers of atoms of each type, C's have a 2-1 abundance ratio, DO's are 3-1, etc. 
Unfortimately, this scheme breaks down as early as A15 . Numbers were assigned in roughly 
the historical order of the study of the lattice. 

• Barrett and Massalski give a table (p. 29) which gives the correspondence between the 
Strukturbericht types and the complexity of the crystal: 



Strukturbericht Designation 


Crystal Type 


A 


Elements 


B 


AB compounds 


C 


AB2 compounds 


D 


Aj^Bj^ compounds 


E,F,G,H-K 


More complex compounds 


L 


Alloys 


O 


Organic compounds 


S 


Silicates 



Strukturbericht Types: 




A T y pe gType CType 
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L'2 (ThH2^ Structure 



LIq (AuCu) Structure 




rCuPt) Structure 





LI2 (Cu j Au) Structure 



L2^ (Heusler) (AlCu^Mn) 
Structure 




L6 q (CuTij) Structure 
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Reference Date: 5 August 1997 
Last Modified: 21 Oct 2004 



The CuAu (Ll^) Structure 




You can now 

• see the structure from several perspectives : 

• examine the structure with an external viewer : or 

• download the coordinates of the atoms in these pictures in X YZ format. 

• visualize the structure (Uses the JMOL Applet) 



• This is a tetragonal distortion of the fee structure. The atoms are at the positions of a face 
centered cubic lattice when c = a. Indium has a related structure. The B19 (AuCd) 
structure reduces to LIq in certain cases, as does the L6 q (CuTi ^) structure. 

• Note that Pearson's Handbook sometimes assigns a Pearson symbol of tP4, rather than tP2. 
Since Pearson's Handbook also says CuTi^ is in the CuAu structure, I'll use my best judgment. 



• Prototype: AuCu 

• Pearson Symbol: ff2 

• Strukturbericht Designation: LIq 

• Space Group: P4/mmm (Cartesian and lattice coordinate listings available) 

• Nvunber: 123 

• Other Compounds with this Structure: TiAl, Nlln 

• Primitive Vectors: 
Aj = y2aX- ViaY 

A2 = VidiX^VidiY 

A3 = cZ 

• Basis Vectors: 
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Bj = 0 (Au) (la) 

B2 = /2Ai + ;/2A2 + /2A3 = ViaX+VicZ (Cu) (Id) 

See these vectors in LaTeX output format 



Go back to the Cubic Close Packed Structurepage. 
Go back to Crystal Lattice Structure page. 
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MECHANISM AND MODELING OF SAW-TOOTH 
STRUCTURE FORMATION IN THE Ll24llo)rWO-PHASE 

SYSTEM' 
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USA and 'Laboraloire d'Etude des Microstructurcs, CNRS-ONERA, B P. 72, 92322. ChfitiUon, Ccdex, 

France 

(Received 9 September 1999: accepted 23 December 1999) 

Abstract— The strain-accommodating mechanism of formation of the saw-tooth microstructure in the 
H2+LI0 two-phase ordered system is proposed. To describe the atomic scale kinetics of ordering and de- 
composition in this system, the master equation, which expliciUy incorporates the effect of the transform- 
ation strain, is formulated. It is used to simulate the predpitation of the Llo phase from the Lh parent 
phase in the Co-Pt alloy. The computer sunulation shows that the decomposition occurs heterogeneously 
on antiphase boundaries of the LU structure if the compoation is near the solubiUty Umit of the LI2 
phase. It eventually produces the saw-tooth microstructure. The decomposition transforms the (OiO) anu- 
phase boundary into a layer of the single-variant tetragonal Llo phase separating the LI2 antiphase 
domains. Later, the strain accommodation transforms this single-variant Llo phase layer into the saw- 
tooth pattern. The simulation results an in excellent qualitative and quantitative agreement with our TEM 
images of Coja sPtii s. © 2(W0 Acta MetaUurgiea Inc. hMished by Eisevier Science Ud. All rights reserved. 

Keywords: Phase Uansformations; Microstructure; Elastic properties; Theory & modeling; Alloys 




1. INTRODUCnON 

TEM observations of the decomposition in Co-Pt 
alloys [1. 2] have shown that the two-phase micro- 
structure in these alloys fonns a very unusual saw- 
tooth-Hke morphology. This morphology has been 
reported in Rcfs [1, 2] for the Cosa.sPtfii.s alloy, 
gradually cooled to 873 K. At 873 K. this alloy is 
in an LI2+LI0 two-phase state with a small volume 
fraction of the Llo phase. As shown in Ref. [2]. the 
cooling initially produces large LI 2 antiphase 
domains separated by thin single-domain layers of 
the Llo phase. These layers ate onpoth as long as 
their width does not exceed abont 3-4 nm. When 
the width exceeds this critical thickness, the layers 
become unstable and suddenly adopt a saw-tooth- 
like shape (see Fig. 1). 

Co-Pt alloys are not the only type of system 
where the saw-tooth structure is observed; it was 
also observed in Cu-Bc alloys by Cheong et aL [3]. 
The structure is formed during the decomposition 
of a hyper-eutcctoid P (b.cc.) phase producing the 
y{B2)+a(f.c.c.) two-phase mixture. The crystallo- 
graphy and kinetics of the decomposition reaction 
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here is much more complex than in the Co38.5Pt«i.5 
alloy. The decomposition in the Cu-Be alloy is ac- 
companied by the B2 — * f.c.c. crystal lattice re- 
arrangement described by the tetragonal Bain 
distortion. The common feature of these two very 
different systems is that the parent phase is ordered 
and that the stress-free transformation strain is 
tetragonal. This circumstance hints that in both 
cases the saw-tooth microstructure is caused by the 
transformation-induced strain and somehow rehited 
to the ordered character of the parent phase. 

The first step towards understanding the mechan- 
ism of the saw-tooth microstructures has been made 
by analyzing the dark field and high resolution elec- 
tron microscope unages of the Co-Pt alloys [2]. It 
has been shown that the LI2 phase has no anti- 
phase boundaries: the antiphase LI2 domains are 
always separated by an Llo phase layer. This is in 
agreement with the prediction that the absence of 
the antiphase boundaries (APBs) in an ordered sys- 
tem is an indication that the ordered phase is 
formed within the two-phase field of the phase dia- 
gram [4]. It has also been shown that the tetragonal 
axis of the LIq layer is always normal to the anti- 
phase vector relating the two Llj domains [1, 2]. 

A detailed mean field study [I, 5] of the APB- 
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enhanced predpitation of the Llo ordered phase 
has shown that the nearest and next-neaicst inter- 
action is sufficient to explain both the formation of 
the Llo layer on the APBs and the relationship 
between the directions of the tetragonal axis and 
antiphase shift vectors. However, this short-range 
interaction model is unable to explain the formation 
of the saw-tooth morphology. The purpose of this 
paper is twofold. We developed the kinetic model 
that incoTporates the strain energy into the thermo- 
dynamics and kinetics of the coherent decompo- 
sition so that the model is applicable irrespective of 
the geometrical complexity of the microstructure 
pattern. Using this model, we reveal the mechanism 
of formation of the saw-tooth structure. It will be 
shown that: 

1. if the alloy composition is in a vicinity of the 
ordered phase solved on the phase diagram, the 
precipitation of the second phase mainly occurs 
on antiphase boundaries; 

2. the saw-tooth morphology develops as a result of 
a relaxation of the long-range coherency strain. 



2. THEORFnCAL MODEL 



The faystal Ij^l^ misfit between the cubic LI2 
and tet^^ettft}mn)phases is described by the tetra- 
gonal stress-free (nmsfonnation strain tensor. The 
Llo phase forms three tetragonal orientation var- 
iants in the fee. structure whose unit cells are 
shown at the top of Fig. 1(a). Each variant gener- 
ates its own stress-free tensor. These tensors arc re- 
lated by the symmetry operations of tiie parent 
f.c.c. phase. The misfit between lattices of the orien- 




Fig. I. Dark field image of the Co3s.sPt6i.s alloys graduaUy 
cooled from 1023 to 873 K at the rate of 40 K, aged for 
IS days at 873 K and water quenched. The dark field 
image was obtained in a zone axis close to (001) with the 
superstructure reflection (210). 



tation variants and parent phase generates the elas- 
tic displacraeots and related ekistic strain. The 
strain energy caused by the phase transformation is 
stored in the strain field. Unlike the chemical free 
energy that depends only on the volume of the 
phases, the strain energy also depends on the shape 
and spatial pattern formed by the Llo phase pre- 
cipitates. Under these conditions, the strain energy 
is affected by the geometry of the microstructure 
pattern and thus controls the two-phase microstruc- 
ture evolution. This coherency strain is a quite com- 
mon phenomenon for practicaQy all multiphase 
coherent systems. Therefore, any realistic theory of 
the microstructure evohition should Uke the trans- 
formation-induced strain into account 

In systems such as alloys with the LI0 + LI2 two- 
phase field, the strain should be explicitly incorpor- 
ated into the thermodynamics and kinetics of the 
phase transformation. One of the most impressive 
examples illustrating the influence of the strain- 
induced interactions is a formation of the chess- 
board-like microstructures during a precipitation of 
a tetragonal phase from the cubic parent phase. 
The chessboard-like microstructures have been 
observed in many systems such as Co-Pt [1, 2) and 
CuAu-Pt The theory of the chessboard struc- 
tures was developed by Le Bouar et td. [7]. The 
computer simulation in [7] has demonstrated that 
the chessboard structure is fonned as a result of the 
strain-induced microstructure evolution in the cubic 
-+ cubic + tetragonal transformation. 

The transformation-induced strain is also a key 
factor in developing the saw-tooth structure. Special 
conditions under which the saw-tooth structure is 
fonned are actually revealed by the TEM obser- 
vations [1]. According to these observations, the 
saw-tooth structure is fonned on the APBs within 
the LI2 ordered i^se. To describe the processes 
developing on the APBs, the kinetic theory should 
be microscopic. There are two microscopic kinetic 
theories that can be used. One of these theories, for- 
mulated in 18, 9], is appUcable to the case of small 
transformation driving force, i.e. when the system is 
not too far from equilibrium (when the supeisatura- 
tion is small): the evolution rate is assumed to be a 
linear function of the thennodynamic driving force. 
In particular, this approach is applicable for the 
coarsening regime. The microscopic theory [8, 9] 
has been extended in [10] by including the strain 
energy. The second approach to this problem is 
based on the kinetic master equation (ME). 
Although the ME approach is more complicated, 
its applicability can be extended to a regime where 
the transformation driving force is high. This is a 
high supersaturation regime of the growth stage. In 
the latter case, a mean field approximation of the 
ME results in a set of coupled kinetic equations 
describing the evolution of the occupation probabil- 
ities and multiparticle correlation functions, the 
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evolution rate a non-linear function of the thernio- 
dynamic driving force. 

2 J, Kinetic equation 

A new and very effective method of using the 
ME for the nearest-neighbor jump diffusion in a 
binary substitutional alloy has been developed in 
Refs (11, 12]. The method uses an approximation 
based on decoupling the averages of products of oc- 
cupation numbers c(r) related to different crystal 
lattice sites [c(r) is a value equal to 1 if the site r is 
occupied by solute atom and 0 otherwise]. This 
decoupling procedure results in the mean field ap- 
proximation, which drastically simplifies the ME 
and makes it a convenient tool for the description 
of the microscopic kinetics far from equilibrium. 

If the mixing energies are assumed to be pairwise, 
the corresponding ME equation is 

-„(r)KrV^^^-«(r)(l 

-„(r'))e^^«^l (1) 
where fi° Xjk^T, i^a is the Boltzmann constant, 

Wr))= X;fK(r-r''Wr'') (2) 

is the average local potential at the site r caused by 
the "chemical** short-range interaction with all 
atoms around the site r, and « ^ mixing 

energy of a pair of atoms in the sites r and r". The 
sum in equation (2) is taken over all sites of the 
lattke except r. In equation (1), i(r-r') are the 
transition rate coeffidents defined as 

f/oe-^^ if ir^r'|=a 

[o otherwise 

where /o is the frequency of "attempts" of a diffu- 
sion jump, IJ is the diffusion tiarrier, a is the near- 
est-neighbor distance and z is the coordination 
number. 

Following Ref. [13], the kinetic equation (2) can 
be presented as 

^ =2 J^Kr-r'Kd -»(r)MrXr) 
dl 

d-nCr')))'^ (4) 
where pit) = (0(r)) + /r^ lnWr)/[l - «(r)]) is in 



fact the chemical potential 6F/dn(r) calculated from 
the mean field firee energy of the non-equilibrium 
state: 

-r'E"(') ln(«(r)) + (l-n(r)) ta(l-Kr)) 

(5) 



2 J, Occupation probabiiities and the l.r,o. 
parameters 

The system that we investigated consists of 
pliases formed by ordering over f.cc. lattice sites 
and decomposition into ordered cubic LI2 (CuaAu 
type) and tetragonal Llo structure (CuAu type) 
phases. The preferred positions of the solute and 
solvent atoms in unit cells of the LI2 phase and in 
three orientation variants of the LIq phases are 
shown in Fig. 2. Occupation probabilities n(t) of 
both ordered phases and their translational and 
orientation variants are described by the same 
equation: 



ii(r) = c + iTi exp(iki r) + ri2 exp(ik2r) 
-I- 1^3 exppkir) 



(6) 



where c is the average composition of the alloy, 
r"*ao(x, y, z) is the f.c.c. lattice site vector, <2o is the 



iooii 




0>) 



'1 *i ' i 



J 



Jnit cells of the(lloj^nd LI2 s tructures and their 
on on the (001) fee plane. (a)' Unit cells of three 
variants of the LIq ordered structure, (b) A unit cell of the 
LI} ordered stn^ure. Black circles indicate the solute B 
atoms, open drcles indicate solvent A atoms. Bold lines 
describe a unit cell Dashed vertical lines show the projec- 
tions of alternating columns of B and A atoms along the 
[001] direction on the (001) plane. Intersection of oolunuis 
with the (001) plane forms the 20 square lattice. Small 
black and open circles in sites of the 2D square lattice 
characterize a type of the columns. 
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f.c.c. lattice parameter, (jc. y, z) are coordinates of 
the f.cc. lattice sites (x, y, z are integer or half-inte- 
gers and j: + >' + z = integer); k, = 27i/jo(100), 
k2=2«/flo(010) and k3-2w/ao(00i) are the wave 
vectors belonging to the same star of the superlat- 
tice vectors {100}; I2, m are the components of 
the Lr.o. parameter. Equation (6) provides a re- 
lation between the microscopic parameters, c, l.r.o. 
parameters »7a, and the microscopic par- 
ameters, which are the occupation probabilities R(r). 

Equation (6) describes the occupation probabil- 
ities as a superposition of three concentration waves 
with the wave vectors, ku ki, ka [9, 14]. The coin- 
ponents nu ^27 ^li of the l.r.o. parameter are ampli- 
tudes of three concentration waves. As follows from 
(6), if = " ';3 = 0, equation (6) describes a disor- 
dered state since all the sites are occupied with the 
same probability c. Otherwise, equation (Q 
describes an ordered phase. 

Depending on the composition c, the interatomic 
interactions in the relevant system are such that the 
global minimum of the free energy F(c, tji, 1/2* li) 
degenerated and falls at either of six points [tfi, fjj* 
= to. 0, 0) = ±;i(l, 0, 0), ±m I. 0), ±»?(0. 0, 
1) or at four points {fju niy '?} = »7(^ 1> 

1). 1(1. 1 i). ^i. I J) and n(h i. I). This degener- 
acy of the global minimimi automatically follows 
from the symmetry of the free energy F{C, riu r\2> ni> 
with respect to the symmetry operations of the f.cc. 
latttee. 

Each of six degenerate points {ij, 0, 0}, being sub- 
stituted in equation (6), describes one of six possible 
Llo orientation and translation variants. The LIq 
structure (see Fig. 2) is formed by alternating cubic 
planes filled by constituent components. The tetra- 
gonal axis of each variant of this structure is paral- 
lel to the wave vector of the wave with the non-zero 
amplitude. Each of four degenerate points {r\, ri,t\}, 
being substituted in equation (6), describes one of 
four translation variants (antiphase domains) of the 
Llj phase. Geometrically these antiphase domains 
are separated by three kinds of APB, each of them 
being characterized by the antiphase shift vectors, 
R, = i[Oll]. R2=:J[101] or R3 = J[nOJ. Because 
the concentration wave representation is valid 
for all variants and all compositions of both phases, 
LI2 and Llo» this approach is well suited for a 
characterization of the LI2 — ► Lh+Llo decompo- 
sition and, in particular, it is valid for a characteriz- 
ation of decomposition in Co-Pt. 

Equation (6) is applicable not only to describe a 
structurally homogeneous phase, but it is also valid 
to describe a structurally and compositionally in- 
hoihogeneous mesoscopic microstructurc with the 
LI 2 and Llo ordering. In this case, the concen- 
tration, c(r), and three l.r.o. parameters, ly^r), ;;2(i'), 
i;3(r), are "smooth" functions of the coordinate r 
rather than constants. The term "smooth" implies 
that these fiinctions appreciably change over the 
distances, considerably exceeding the crystal lattice 



parameter. The atomic-scale structure of inhomo- 
geneous mesoscopic microstructure is still described 
by equation (6) uang c = c(r). i?i=i?i(r). »f2'''?2(r). 
fy3 = fj3{r) for each site r. The mesoscale Lr.o. par- 
ameter profiles ripij) and the composition profile 
c(r) can also be, in turn, expressed throu^ the 
microscopic parameters, occupation probabilities 
/i(r), by using the equations 



^r) = i£«(r + a;) C7b) 

where the summation is carried out over the tetra- 
hedron of the nearest neighbors in the f.c.c. lattice, 
the site r being its origin. In (7), «i=*0 and three 
vectors, 82, as and 84, detennine the three nearest- 
neighbor sites of the f.cc. lattice tetrahedron with 
the origin at r The vectors k^ arc three superlattice 
vectors defined in equation (6). Therefore, equations 
(6) and (7) provide a relation between the meso- 
scopic and microscopic scales of the inhomogeneous 
microstructure. 

23. Strain energy of the structurally inhomogeneous 
coherent system 

No realistic theoretical characterization of a 
coherent structural transformation is possible unless 
the transformation-induced elastic strain is taken 
into account. The strain is generated by the crystal 
lattice mismatch between the low symmetry phase 
and the high symmetry matrix, as well as between 
different orientation variants of the low symmetry 
phase. The interaction between finite elements of 
the precipitate phase caused by the strain has a 
dipole-like long-range asymptotic behavior, which 
drastically distinguishes it from the short-range 
"chemical" mixing energies. Because this is the 
long-distance part of the strain-induced interaction 
potentials that affects the microstructure evolution, 
we can use the continuum elasticity for its descrip- 
tion. In fact, the potential calculated by using the 
continuum elasticity is a long-range asymptote of 
the atomistic potential. 

It is important that using the continuum approxi- 
mation does not introduce an inaccuracy in the 
description of the morphology evolution. Indeed, 
the atomistic corrections to the continuum elasticity 
potentials caused by the discrete nature of a crystal, 
by definition, are substantial at distances commen- 
surate with interatomic distances and thus are short 
ranged. Therefore, these atomistic corrections can 
be incorporated into the short-range "chemical" in- 
teraction potentials. As a result, the short-range 
energies of mixing, entering (S) become tenor- 
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malizcd effective mixing energies. They include the 
short-range atomistic corrections to the continuum 
elasticity potential. This procedure, in effect, elimin- 
ates all inaccuracy in the use of the continuum elas- 
ticity instead of crystal lattice statics because in this 
approach all corrections to the continuum elasticity 
prove to be incorporated in the mixing energies W{, 
which arc fitting parameters anyway. They are 
determined by fitting the calculated phase diagram 
to the observed one. 

The remaining strain interaction potentials are 
described in terms of the continuum elastidty. 
These potentials are calculated below by using the 
Khachaturyan-Shatalov (KS) Phase Fidd 
Microdastidty theory of structurally inhomo- 
geneous systems [9. 1^. The total free energy of the 
system is a sijm of the "chemical** free energy (5) (it 
also iochides the short-range corrections to the con- 
tinuum elastidty) and the elastic energy, i.e. 

The KS theory provides the analytical solution of 
the elastidty problem for a body that is macro- 
scopically homogeneous but structurally inhomo- 
geneous in the mesoscopic scale. This is a situation 
where the typical inhomogendty length is much 
smaller than the size of the body. In particular, the 
KS Phase Field Microdastidty theory gives the 
closed analytical equation for the strain energy 
function of the continuous concentration and l.r.o. 
profiles (phase fields), characterizing an arbitrary 
coherent microstructure. Using this analytical sol- 
ution makes it unnecessary to calculate a numerical 
solution of the elasticity equation for each time 
increment of the microstructure evolution. This 
drastically shortens the numerical calculations and 
allows one to address the problems which are 
intractable or very time consuming if the finite d- 
cment methods are used. The KS theory was in- 
itially applied to particular step-like phase fields 
describing structurally homogeneous coherent in- 
clusions with sharp boundaries. The theory has 
been developed further by Chen et al. [16], Wang et 
ai [17] and Lc Bouar et al. [7], who applied this 
formalism to the cases of arbitrary continuous 
phase fields, the concentration and l.r.o. order par- 
ameter fields. 

In the case where the strain is predominantly 
caused by the l.r.o. parameter heterogendties, the 
stress-firee strain determining the crystal lattice 
misfit at a point r can be presented as a function of 
72 and 173: = fh. 'Js). The stress-free 
strain is invariant with respect to rigid-body trans- 
lations which produce the transf<mnations of the 
t.r.o. parameters: 

(fl}* »?2» f3)"*(-'?l» '?2» - 73). ^2* ^i) (8) 



Expanding ej/fji, rji, ri^) in the Taylor scries up to 
the first non-vanishing terms and eliminating the 
expansion terms that are not invariant with respect 
to the transformations (8) leaves us with an 
equation which is invariant with respect to any 
rigid-body translations. The general form of the 
first non-vanishing term of such an expansion is: 

where the tensor ^{p) is the strain transforming 
the crystal lattice parameten of the f.cc. phase into 
the parameters of the equilibrium pth variant of the 
low symmetry phase, and m is the total number of 

variants. 

In the case of an LI 2 Llo transformation, the 
number of possible orientation variants is m » 3 
and the stress-free strain (9) can be presented as: 

„2/c3 0 0 \ „2/ei 0 0\ 

''o\o 0 e,/ "^XO 0 e,y 

0 0\ 

+ ^4 0 a, 0 (10) 

^0 \o 0 ej / 

where = (c - flo)/flo» s\={a — oo)/ao, «d 
{a, c) are the ciystal lattice parameters of the cubic 
f.cc. and the tetragonal Llo phases, and tio is the 
equilibriim l.r.o. parameter of the LIq phase. 

Following Ref. [16], we can express the strain 
energy as a functional of the Lr.o. fields: 

where {»jj(r)}fc is the Fourier transform of i7j(r), 
(i}|(r)}][ is the complex conjugate of {rii(r)}k and 
B„(^) = C^i^{p)JS(9)-et^^py^(€^^^^ 
e «k/Jt is a unit vector in the redprocal space and e« 
its ith component, djki is the elastic modufi tensor, 
<^^p)=C(fia^ip), and Sli^{e) is a Green function 
tensor which is inverse to the tensor Q^* = Cjki/ekeh 
The integration has to be carried out over the first 
Brillouin zone, where a volume of {In^jV about 
k=0 is excluded from the integration. When the 
vohime V is large, this defines the prindpal value of 
the integral. 

To make MB equation (1) suitable for character- 
ization of the transformation producing the coher- 
ency strain, we have to modify the local potential 
(0(r)) by adding to it the contribution (0(r)ei) gen- 
erated by the strain. The latter transforms equation 
(2) into: 

{<m = jyi^ - )«(r" ) + {«»ei(r)) (12) 
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where (0(r)ei> is the variational derivative of (11) in 
R(r). It can be rewritten as: 



{Oei(r)> = 



where the vectors are the three supcrlatticc vec- 
tors defined in equation (6) and the vectors ay are 
defined in equation (7). 

For numerical calculations it is convenient to 
express the integrals over k, as sums over all N vec- 
tors k within the first Brillouin zone are permitted 
by the periodical boundary conditions. This gives: 



«»d(r)): 



, p=3«=3>-4 

4e:*«Q)«'^«»»«^'^' 

where v s is the vohime per f.cc. lattice site, 
N is the total number of the f.cc. lattice sites, the 
prime sign means that the term k«0 is excluded 
from the summation, 

l»lJ(r))»=E''?«"P(*')- 



This transition from integration to summation in 
equation (IS) and (16) has the macroscopic accu- 
racy of the order of l/N. 

The driving force (14) is a functional of the I.r.o. 
parameters rjp(ty Using equation (7) in (14). we 
express the ME (1) in terms of the function /i(r). 
The obtained kinetic equation with respect to 
n(r, 0 nonlinear finite-difference differential 

equation whose solution fully describes the micro- 
scopic and mesoscopic microstructure. Computer 
simulation of the microstructure evolution is 
reduced to the numerical solution of this equation. 



3. HETEROGENEOUS DECOMPOSITION NEAR 
ANTIPHASE BOUNDARIES 

To insure that the free energy (5) provides the 
phase diagram with the Lh+Llo two-phase field, a 
proper interaction model has to be chosen. The sim- 
plest short-range interaction model which results in 
the Lh+Llo two-phase conccnUation range is the 
nearest and next-nearest interaction model, with the 
positive nearest-neighbor mixing energy FTj and 
negative next-nearest energy Wz. In this work we 



consider the case Wj/fVx^-OX > 0- At the 
temperature T = XMSWi/kji, the calculated 
Llj+Llo two-phase field is within the range 0.42 < 
c < 0.43 at this temperature. This range is of the 
same order as the observed LI2+LI0 two-phase 
range in the prototype Co-4H system (approxi- 
mately 0.395 < cct, < 0.405 at 963 K, Ref. [IJ). 

S.L Short-range interaction 

The formation of the saw-tooth structures is 
observed on the APBs of the LI2 phase, which were 
almost perpendicular to a cubic plane. Below wc 
consider the APB in the LI 2 lattice contained in the 
(010) plane. In this case, the observed mesoscopic 
saw-tooth microstructures are homogeneous along 
the [001] axis and the saw-tooth pattern is located 
in the (001) plane. The homogeneity of the hetero- 
geneous microstructure along the [001] axis actually 
means that there is no diffusion along this axis. 
Under this condition, the 3D problem of the micro- 
structure is reduced to the 2D problem. The 
assumption of homogeneity along the [001] direc- 
tion drastically simplifies the simulation model: it 
implies that concentration c and I.r.o. parameters 
7i* ill ^ equations of the 3D model do not 
depend on the coordinate z along the [001] axis, i.e. 
c(r) = c(jc, y) and iip=??p(x, >'), where x and.;* now 
designate the coordinates of columns of atoms in 
the f.c.c. lattice parallel to the [001] direction. The 
. intersection of these columns with the (001) plane 
generates the square lattice whose parameter is <2q/2 
(see Fig. 2), Therefore, the 2D position of the 
atomic columns is r =: ^ a^fyc, y), where the coordi- 
nates x, y of the columns are integers. 

Counting all nearest and next-neatest interactions 
within a column and between the neighboring col- 
umns g^ves us the following mixing energies 
between columns: 



^Fchcm(r-r') = 



2}V2 if|r-r'!=0 

2JVi if |r-T'|=ao/2 

Wi if |r-r'|=flo/V2 (16) 

W2 if|r-rM=flo 

0 otherwise 



where r' = J ao(x', yO- As long as the microstruc- 
ture homogeneity along the [001] axis is maintained, 
this 2D model is mathematically equivalent to a 3D 
model with mixing energies, given by equation (16). 
Because the LI 2, Llo and the f.cc. solid solution 
are compatible with this homogeneity condition, the 
2D model correctly reproduces the equilibrium 
structures, the mean field kinetics, and thermodyn- 
amics of the 3D model. 



i.2. Strain energy 

To establish an origin of saw-tooth microstruc- 
tures and distinguish the role of the different factors 
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in its f6rmation, we exclude the efifect of the elastic 
moduli anisotropy. As a matter of fact, the elastic 
anisotropy is not a factor which considerably com- 
plicates the computations. However, as it has been 
shown m our previous paper {7J, even the clastically 
isotropic model predicts the chess-board structure 
in Co-Pt and other alloys, reproducing with great 
accuracy the finest features of their microstructure 
morphology. It is shown below that the clastically 
isotropic model also conectly predicts the for- 
mation of the saw-tooth morphology and its mor- 
phologic features with excellent quantitative 
agreement with the observations. This actually 
proves that it is the crystallography of the crystal 
lattice rearrangement (symmetry and the relation 
between components of the misfit tensor) rather 
than the elastic moduli anisotropy that is the major 
factor detennining the mechanism of both the saw- 
tooth and chess-board structures. 

Homogeneity along the [001] axis simplifies 
equations (ll)» (13) and (14) for the strain-induced 
driving force. This homogeneity results in vanishing 
of in all equations involving the strain energy. 
Therefore, we used k'^ik^, ky, 0) and e^(ex, Cy, 0) 
in all equations involving the strain energy because 
all integration (summation) over is lifted. With 
this simplification^ the assumption of the elastic iso- 
tropy, and the stress^free strain (10), the functions 
BpqWk) defined in equation (U) are: 

iJ|l(«x,«y>0) 

= 7tM-j(«i - *)^« - + *^))' + (1 - 
(1 - 



533{e„ey,0) = 



(l+2vdb + rfo) 



B|2(Cx.«/.0) 



the strain-induced driving force for arbitrary l.r.o. 
parameter profiles. 

3 J. Master equation 

As has been discussed in Section 3.1» there is no 
diffusion along the [001] direction if the decompo- 
sition maintains homogeneity along this direction. 
Then substitutional diffusion develops between the 
[001] colunms in the f.c.c. lattice and can be 
described in terms of the 2D diffusion over the 
primitive square lattice with the oofl parameter. 
Master equation (2) for the diffusion in the f.c.c. 
lattice is equally accurate for the 2D model of the 
diffusion in the square lattice if the elementary dif- 
fusion events are the nearest-neighbor jumps and 
the following probafaiiities for jumps between col- 
umns are assumed: 



= - *>'^ • ^ + *(1 + *X1 + V)) 



= ^j^^((4-l-v)(l-*)^J + *(*+!) (17) 
x(l+v)) 

where e is a unit vector, v and ft are the Poisson 
ratio and the shear modulus, respectively, and 
db'-fii/es is the stress-fi-ee strain parameter. The 
functions ^22(^jc, 0) and B^^iext Cy, 0) are 
obtained from Bn(«x. ^y. 0) and Biiie^,, Cy, 0) in 
(17) by the simultaneous and XT^y cyclic per- 
mutation. Using equation (17) in equation (14) 
allows us to formulate the closed analytical form of 



Kr-r'). 



If if |r-rM=flo/2 

/ if lr-r'|=flft/V5 

-12/ ifr-r' = 0 

0 otherwise 



(18) 



where r = ^ a^{x, >') and r' = ^ ao(^', y') and x, y, 
and x\ y' are integer coordinates, /= 
/o &Lp{-U/(kBT)) is the jump frequency between 
the nearest-neighbor sites of the 3D f.cc. lattice. 

It is convenient to solve equation (1) using the 
reduced variables and parameters. The reduced time 
is r*«/r and the reduced temperature is = 
khT/Wi. The strain energy contribution is deter- 
mined by the dimensionless parameter 



(19) 



where v is the volmne per f.cc, lattice site and x 
determines a balance between the strain and 
^'chemical" energy of the alloy. 

3.4. Input parameters for Co-Pt alloys 

To apply our simulation model to the Co-Pt 
alloys, we have to estimate the values of the input 
parameters, using for this the data obtained in the 
measurements tmrdated to the data obtained from 
the geometry of the microstructure patterns. A 
choice of the mixing enei:gies should provide the 
ordering that produces the LI 2 and Llo phases and 
their two-phase equilibrium. The two-phase concen- 
tration range should be of the order of 1 at.% and 
should be situated near Cco'^0.4. These are the esti- 
mated values that follow from the TEM obser- 
vations of the Co-Pt alloy at 963 K [2]. At 
Wz/Wi = -0.2 and T* = k^T/W^ = 1.125 at 
> 0, our model gives the LI2 and Lip two- 
phase range which is close to these estimates. It 
does automatically produce the LI2 and LIq phases 
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that coexist in equUibrium within the 0.4215 < c < 
0.4320 compositioD range. 

The most important data for the microstructure 
prediction are the crystal lattice parameters of the 
phases at different temperatures and compositions. 
These data determine the strain energy contribution 
and all specific features of the two-phase micro- 
structure patterns. The crystal lattice parameter 
data allow us to calculate the strain component 
entering x equation (19) and the parameter do 
determining the anisotropy of the function Bpg(€) 
in (17) characterinng the strain-induccd interaction 
between elements of the microstructure. 

The crystal lattice parameters that we need to 
know are the parameter of the fee. disordered 
phase, Ofxx.* ^ parameter of the cubic ordered 
phase, a(Ll^, and the parameters of the tetragonal 
Llo phase, Ct and at* within the LI2+LI0 compo- 
sition range at the decompodtion temperature (at 
cco ~ 0.4 and T = 873 K). Unfortunately, these 
data are not available for those compositions and 
temperatures where the observed saw-tooth struc- 
ture has been formed. Below we assume that 
ii^-. 61/63 =0.2 and £3 =-0.02, The following con- 
sideration indicates that these values provide a 
reasonable estimate. 

Indeed, according to equation (10), the tetragonal 
ratio cjot of the Llo pbase within the two-phase 
region, expressed through the constants et<»^3 
azKl 63, is 



With dS,=0.2 and 63 = -0.02, equation (20) gives ct/ 
flt= 0.984. This value is consistent with the exper- 
imentally observed tetragonalities of the Llo pbasc: 
for example, the observed tetragonality of the 
C047.6Pl5i4 aUoy is cja,=0,m [18]. The tetragon- 
ality ct/flt= 0.984 following from the assumed input 
parameters is closer to unity than the observed tet- 
ragonality. c/flf- 0.973. However, this deviation is 
not an indication of a disagreement with obser- 
vation results: the composition considered in our 
model, cco^OAy is farther off from the stoichi- 
ometry, cco'-O.S, than the composition, 
ccb = 0.476, of the alloy whose parameter is 
reported in [18], and thus the tetragonahty ratio of 
our aUoy is expected to be closer to unity. 

The temperature of our TEM study of the saw- 
tooth structure in Co-Pt alloys was T - 873 K. 
Our choice, r-^B^/'^i = 1125, providing the 
proper concentration range of the LI2+LI0 field 
corresponds to this temperature if Wxlk^~ll^ K. 
This value is used in our simulation. 

The estimate fKi//rB='776 K is obtained by balan- 
cing the mean field free energies of the two ordered 
phases. Since the mean field theory is more accurate 
for ordered phases, this estimate is more reliable 
than the one obtained from the f.cc. — > LI© order- 



disorder transition temperature Tq. For the tran- 
sition temperature To, the mean field theory gives 
ro= 1.375 W\lk^ at c=0.5. It results in FKi/^b = 
798 K if the experimental value for the transition 
temperature, To^lO^S K. is used. The latter esti- 
mate is less reliable than the pre^ous one because 
the area near the order-disorder transition point is 
an area of the phase diagram where the mean fidd 
theory is least accurate (it usually gives the order- 
disorder temperature as 20-30% higher than that 
calculated by more accurate methods). However, it 
is still in line with the estimate W\lk^ = 776 K. The 
equilibrium value of the l.r.o, parameter of the Llg 
phase within the LI2+LI0 region following from 
our model at F^,/ikB=»776 K is ?/o= 0.2685. 

The estimate for the shear moduhjs of the Co-Pt 
alloy gives /i 70 GPa (interpolation between the 
value ^=81 OPa for pure Co and ^=62.5 GPa for 
pure Pt). Using the parameters /4 ~ 70 GPa, 
iro-0.2685, e3=-0.02, Wxlk^^m K, v = aJ/4. 
oq" 3.8 A in (20) gives the estimate for the input 
parameter 

70 X 10^0.02)^ (3.8 X lO-^^)^ 

' 776* 1.38 X 10-"* (0.2685)* 4 

Therefore, onr input parameters are: 

<4 = 0.2. x = 7. r* = 1.125 (21) 



4. SIMULATION BESULT5 

To find the microstructure evolution as a funo> 
tion of time the non-linear differential equation 
(1) for n{x, y, O was solved at each point r*=(x, 
y) by using the Buler method. A small random field 
was added to avoid trapping in local minima. The 
initial occupation probabilities, n(x, y» /'-O). were 
chosen to describe two unrelaxed LI 2 antiphase 
domains with homogeneous composition c« 0.4235. 
The APB was initially situated in the (010) plane 
and its antiphase vector was ^[011]. A simulation 
box contained 128 x 128 sites (the [001] columns in 
the fee. lattice) that form the square lattice with 
the crystal lattice parameter flo/2. In order to use 
periodic boundary conditions in all directions, we 
have introduced two antiphase boundaries assod- 
ated with opposite translation vectors. 

The dependence of n^c, y, O vs time /" obtained 
by the numerical solution of the kinetic equation 
(I) simulates the subsequent stages of the mioro- 
structuie development in the heterogeneous LI2 -+ 
LI0+LI2 decomposition developing on the anti- 
phase boundary. The simulation results are shown 
in Fig, 3. In Fig. 3, the two-phase morphology is 
characterized by the coarse-grain concentration pro- 
file c(x, y). This profile is visualized by shades of 
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gray. The white contrast describes the Llo phase 
[greater value of c(x, y)l For a corapaiison, we 
also presented in Fig. 3(d) the micrographs of our 
TEM observation of Co38.5Pt6i.5 annealed at 873 K. 
Figure 3(a) shows that the APB located along 



the (010) plane remains flat, but at /* = 100 they 
arc consumed by the Llo precipitate. The precipi- 
tate is shown by a white contrast. At this point, our 
simulation results are still similar to those obtained 
in Refs [1, S]. where no elastic effects are taken into 
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Fig 3 Comparison between simulation results (aHc) and a TEM image of a saw-tooth morphology. 
The concentration maps (aHc) correspond to the reduced times: (a) r*- 100; (b) r-200; (c) ( =500, 
The dark field image (d) of CojasP^MS was obtained in a zone axis dose to (OOl) with the superstruc- 
ture reSecUon (120). The alloy sample was graduaUy cooled from 1023 to 873 K at the rate of 40 K, 
aged for 15 days at 873 K and water quenched. 
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account The Llo layer consuming the APB bound- 
ary has the (010) habit plane coinciding with the 
APB plane orientation. The Llo phase layer has the 
tetragonal axis parallel to the [100] direction. 
Therefore, the tetragonal axis is perpendicular to 
the ^fOUl out-of-phasc vector of the APB. This 
result is in complete agreement with the experimen- 
tal observations in the Co-Pt system [1, 2]. 

At r*-300 [Fig, 3(b)], the situation changes. 
The precipitate with the (010) habit starts to 
evolve adopting a wavy shape, which indicates 
that the Llo phase plate tends to deviate its habit 
plane from the (OlO)f.c.c.. 

Finally, at f'^500 [Fig, 3(c)], the Llo phase 
forms a saw-tooth morphology. Comparison of 
Figs 3(c) and (d) illustrates the excellent agreement 
between our simulation and experimental results for 
the C6-Pt system. 

It is interesting to note the diffuse contrast that is 
mostly localized around the edges <^ the saw-tooth. 
This occufs because the elastic strain around the 
edges is especially high and deviations of concen- 
tration and Ir.o. parameters from their equilibrium 
values reduce the elastic strain. This phenomenon is 
illustrated in Fig, 4 where an l.r.o. parameter profile 
in a [010] direction is presented. This profile inter- 
sects a tip of the saw-tooth morphology obtained 
when r'' = 500. In the bulk of the Lh domains, the 
three components of the l.r.o. parameter are not 
strictly equal (in absolute value) as they should be 
in the absence of the strain field generating the saw- 
tooth. Furthermore, the value of the second com- 
ponent of the Llo structure is lower than the stress- 
free vahie i|-= 0.2628. 

5. DISCUSSION 

Figure 3 shows that our simulation results are in 
excellent agreement with the TBM observations [1, 
2] of the Co38.s-Pt6i.5 alloy. The sknulation indi- 



cates that the nucleation of the LIq phase occurs by 
the heterogeneous mechanism. The tetragonal Llo 
phase forms at APBs as a single-domain layer that 
gtaduaUy consumes the APB [Fig. 3(a)]. The tetra- 
gonal axis of this layer is perpendicular to the anti- 
phase shift vector of the APB. The antiphase shift 
ao[0| ^] for the (010) APB resulted in the formation 
of the domain of the Llo phase whose tetragonality 
axis is parallel to the [100]. Such domains will be 
called x-domains. The (010) plate of the LIq phase 
[white strip on Fig. 3(a)] is not stable. It cvoWes 
into a saw-tooth structure formed by alternating 
plates of x-donmins with dififerent habit planes [Figs 
3(b) and (c)]. This is also in an agreement with the 
TEM observations [I, 2], 

It should be emphasized that practically no con- 
straints were imposed on the system to choose the 
optimal transformation path from its initial state of 
two LI 2 antiphase domains separated by sharp 
unrclaxed APB. The entire evolution path illus- 
trated by Figs 3(a)-(c) is fully determined by the 
kinetic equation (1). 

The simulated microstnicture shows not only 
qualitative but also excellent quantitative agreement 
with the TEM observations of the CoM^5-Pt6i.5 
alloy. In the emulated picture [Fig. 3(c)], the angle 
between the normals to the tooth surfaces (between 
the habit planes of the x-domain plates forming a 
"tooth" of the structure) and their tetragonal axis 
directions is about 35°. This angle coincides with 
the one observed in the TEM pictures [Fig. 3(d)]. 

In the simulation, the microstnicture evolution is 
practically stopped when the structure shown in 
Fig. 3(c) forms. The typical size of a "tooth** 
obtained at this stage is about 10 nm. A similar size 
is observed in the TEM picture after 15 days of 
aging of the Co38.r-Pt6i.s alloy at 873 K [Fig. 3(d)], 
Therefore, tte scale of the amulated saw-tooth 
microstnicture is m agreement with the TEM obser- 
vations (compare Figs 3(c) and 3(d)]. 
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Fifi, 4. (B) l.r.o. profile along [010] intersecting a tip of the saw-tooth morphology whai t -500. The 
Lb domains of the type r;(l 111 and »?[lH] are separated by a zone of Llo structure of Ac type ijoIlOOJ. 
The path corresponding to the l.r.o. profile (a) is sketched by an arrow in (b). 
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The most challenging part of the theoretical 
characterization of the system is the accurate deter- 
mination of the transformation kinetics in real time. 
The real time t is related to the reduced time C of 
our simulation by equation 

t^r/f=r/{fo exp(-i;/(A:Br)))=» 
i//o cxp(i//(*Br))r (22) 

The coefficient relating real and reduced time in 
(22) is very sensitive to the activation energy of dif- 
fusion U and the aging temperature T. Besides, the 
accurate data on the diffusion coefficients at the 
composition of the simulated Co-Pt alloy are not 
available. To estimate the value of V, we have used 
the activation energy V «= 279.3 kJ/mol for diffu- 
sion of Pi in Co [19]. The typical number of 
"attempts" /o for a djffusional jump is of the order 
of a typical frequency of the crystal lattice vi- 
bration. /o'^ 10'^ s~*. Using these values in (22) for 
/• = 500 at r 873 K gives /-^^SxlO^ss 
30 days. The TEM picture of the saw-tooth struc- 
ture shown in Fig, 3(d) was obtained by gradual 
cooling of C03«.5-Pt<5i.5 from 1023 to 873 K at the 
rate of 40 K, followed by aging for 15 days at 
873 K. Given the fact that (i) the alloy has already 
been partially transformed during gradual cooling 
prior to aging at 873 K and (ii) the input data for 
this estimate are very approximate, the estimated 
structure formation time, 30 days, is in reasonable 
agreement with the 15 days of aging time in the real 
system. In cases like this, an agreement that is bet- 
ter than an order of magnitude can hardly be 
expected. 

All input parameters that are used in the simu- 
lation can be found from independent measure- 
ments of the crystal lattice parameters and the 
phase diagram, i.e. from the data unrelated to the 
morphology of the microstructure. The simulated 
structures are sensitive to the crystal lattice misfit 
data» which determine the crystallographic par- 
ameter do = €i/ey and reflects the geometry of the 
crystal lattice rearrangement in the phase trans- 
formation. The crystal lattice misfit determines the 
magnitude and spatial anisotropy of this strain. The 
sensitivity of the microstructure to the crystallo- 
graphic misfit data follows from the fact that the 
long-range strain is a major factor determining the 
morphology of the mesoscopic microstructure. 

The geometry of the microstnicture is much less 
sensitive to the parameter x weighting the contri- 
bution of the strain energy with respect to the 
"chemical" energy. The simulation is also not very 
sensitive to the choice of the parameters W{r)l 
(tcBT), as long as this choice produces the LI 2 and 
Llo ordering and the correct two-phase range at the 
simulation temperature. 

As is known, the mean field approximation 
underestimates the contribution of the atomic corre- 
lations to the alloy thermodynamics, especially for 



the disordered state. Although the simulated micro- 
structures in this study are obtained by the solution 
of the mean field version of the ME, the agreement 
between the simulated microstructures and the 
TEM observations is excellent. The latter indicates 
that the accuracy in the strain energy treatment 
here is much more important than an underestima- 
tion of atomic correlation. It also shows that the 
KS theory is well suited to treat the strain energy 
contribution. 

The second reason why the mean field ME proves 
to be sufficient to describe this system is that inac- 
curacy of the mean field approximation is partially 
compensated by a choice of the "chemical" energies 
Wu Wi as fitting parameters to the observed phase 
diagram. Besides, as is well known, the higher the 
l.r.o.. the more accurate becomes the mean field ap- 
proximation (it becomes exact at 0 K where the 
l.r.o. reaches the maximum possible values). The 
latter is also a factor since we consider the processes 
where all phases are always m the ordered state 
where the higher accuracy of the mean field theory 
should be also expected. 

The simulation results shown in Figs 3(aHc) 
demonstrate details of the saw-tooth structure for- 
mation. The formation mechanism can be qualitat- 
ively explained by the combined effect of the APB 
and transformation strain. The APBs in the LI2 
phase are high-energy defects. Thus, the free energy 
decrease due to the transformation of APBs into 
the Llo P^ase is greater than it would be if the Llo 
phase was formed in the bulk of the LI 2 parent 
phase. In this situation the Llo phase forms by the 
het^ogeneous nucleation and the APBs are the 
nucleation sites. The Llo phase formed on the 
APBs consumes than, replacing the APBs by a 
single-domain layer. The initial orienution of this 
newly formed Llo inherits the orien- 

tation of the APB, which is the (010) plane in our 
case. This situation is reproduced by the computer 
simulation [Fig. 3(a)]. However, the (010) orien- 
tation is not a strain energy minimizer and thus is 
not stable. The strain energy minimum is achieved 
if the normal ei to a plate of the Llo phase 
"rotates" from its initial direction ejoioj normal to 
(010) to either e,-(cos 35", +8in 35^ 0) or 
C2=(cos 35**, -sin 35"*, 0) directions. The LI© plates 
with both types of habit planes, normal to ei and 
€2, have the same strain energy. 

However, this habit plane rotation occurs only if 
the decrease in the strain energy outweighs the 
increase in the interfadal energy due to a deviation 
of the interface from the minimum interfacial 
energy orientation [9|. It does not occur for thin 
plates since the interfacial energy of thin plates 
plays the dominant role. Thus, the plate is parallel 
to the plane of the minimum interfadal energy. In 
our case, the minimum interfacial energy plane is 
the (010) plane parallel to the APB. When the 
thickness D of the plate exceeds the critical value 
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Da during the growth process, the bulk strain 
energy becomes more important than the interfacial 
energy and the plate starts to rotate to either of 
two orientations dictated by the bulk energy mini- 
mum. A balance between strain and interfadal ener- 
gies of the plate [9] results in the following criterion 
for the critical thickness: 



where AB = [5ii(ei) - -Bii(€[oioi)l < 0 and 
Ay = [y(ei) - Kcpio))) > 0. y(e) is the LIq/LU inter- 
facial energy of the interface normal to the unit vec- 
tor c. If the thickness O of the LIq interfacial layer 
is less than Den the layer just replaces the LI2 anti- 
phase domains and has the (010) habit. During the 
growth, when the thickness D exceeds Den these 
plates start to deform forming the sawtooth struc- 
ture. The teeth of the saw-tooth suiicture is com- 
prised of the plates of the Llo phase, which 
determine a shape of teeth. Points where the nearest 
plates join arc the tooth tips. This transfonnation 
process is illustrated by Figs 3{aHc). The TEM ob- 
servations fully agree with these conclusions. 
According to observation [1, 2], the smooth wetting 
layers of the Llo phase separating the LI 5 phase 
antiphase domains transform into the saw-tooth 
structure when the thickness of the layers exceeds 
3-4 nm. 

In fact, the qualitative consideration presented 
above explains only why the Llo phase plates rotate 
so that their habit is normal to either ej or vec- 
tors when the thickness of the plates exceeds the 
critical value. This is not suffideot to explam why 
these plates rotate to self-assemble themselves in the 
continuous saw-tooth zig-zag pattern. The latter 
can be understood if we take into account that the 
rototions of the elements of a single (010) Llo 
domain plate [Fig. 3(a)] into a set of the LIq plates 
with the different habiU should occur by maintain- 
ing the continuity of this domain. This continuity 
of the domain is required to prevent a formation of 
new LI0/LI2 boundaries on the edges of each 
plates, since the formation of edges would raise 
both interfacial and strain energies. It is an attach- 
ment of these plate, edges that produces the continu- 
ity of the Llo domain, which results in the 
formation of the saw-tooth structure. 

Eventual freezing of the saw-tooth structure is a 
result of extreme difficulties in the coarsening after 
the growth stage is completed. Indeed, the geometry 
of the saw-tooth structure is sudi that its coarsen- 
ing, which maintains the Llo phase continuity and 
orientations of the Llo plates forming the zig-zag 
pattern, would require a collective microstructure 
reconstruction. This process would involve simul- 
taneous and synchronic evolution of several neigh- 
boring plates of different orientations, i.e. the 
process that is much slower than the conventional 



Ostwald ripening of separate particles. Freezing of 
the microstructure evolution observed in our simu- 
lation after the time period 500 is consistent 
with this conclusion. 

Finally, it is very interesting that our simulation 
indicates that the concentration and order par- 
ameters deviate from the equilibrium stress-free 
values to tower the misfit between the Llo layer and 
the LI2 matrix (Fig. 4). In fact, the predse vahies 
of the l.r.o. parameter and the concentration under 
the elastically constrained condition depend on the 
microstructure, and cannot be inferred in advance. 
This effect cannot be predicted by using the sharp 
boundary approximation that usually neglects the 
freedom for the system to relax the l.r.o. parameter 
and concentration profiles from their equilibrium 
values to reduce interfacial and the strain energies. 
An advantage of the field methods is that they do 
not impose any a priori constraint on the way the 
system balances the different contributions to the 
total free raergy. 



& SUMMARY 

This work presents the first microscopic kinetic 
model which describes the major structural charac- 
teristics of the decomposition in the LI2+LI0 stab- 
ility region of the f.cc. binary alloy. The proposed 
theory and the computer simulation explain the 
mechanism of the saw-tooth structtu« formation in 
the Co-Pt generic system and provide insight into 
the conditions where this structure can be expected. 
The saw-tooth structure is formed within the 
Llo + LI 2 area near the solubility line of the Llj 
phase due to the transformation strain accommo- 
dation. 

The model predicts that the LI2 ~* Llo trans- 
formation starts from a formation of Llo phase 
layers wetting the APBs. In agreement with the 
TEM observations, the tetragonal axis of the Llo 
structure is normal to the antiphase vector. When 
the thickness of a layer reaches about 3-4 nm, our 
model predicts the formation of the saw-tooth 
microstructures. The saw-tooth structure does not 
alter the global orientation of the APB but drasti- 
cally reduces the strain energy by a local reorienta- 
tion of the Llo phase layers along the elastically 
soft directions. 

The simulation study demonstrates the ability of 
the model to provide excellent quantitative agree- 
ment with the TEM observations. 
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